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To describe the small-scale intermittency of turbulence, a self-similarity is assumed for the prob¬ 
ability density function of a logarithm of the rate of energy dissipation smoothed over a length scale 
among those in the inertial range. The result is an extension of Kolmogorov’s classical theory in 
1941, i.e., a one-parameter framework where the logarithm obeys some stable distribution. Scaling 
laws are obtained for the dissipation rate and for the two-point velocity difference. They are consis¬ 
tent with theoretical constraints and with the observed scaling laws. Also discussed is the physics 
that determines the value of the parameter. 

PACS numbers: 47.27.Ak 


I. INTRODUCTION 

Consider a fully developed state of three-dimensional 
homogeneous and isotropic turbulence that occurs in an 
incompressible fluid. The kinetic energy is supplied at a 
large length scale and is dissipated into heat by viscosity 
v at around a small scale of the Kolmogorov length 
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Between these two, there is an inertial range of length 
scales, where the kinetic energy is transferred on average 
to the smaller and smaller scales without essentially any 
dissipation. This range has no characteristic scales about 
its lengths r and about random variables defined at each 
of the lengths r. For the spatial average (•) of such a vari¬ 
able X r , we expect the scale invariance {X ir ) = C- y (X r ) 
with C 7 > 0 that is a function of 7 > 0. The solution 
is a power law (X r ) oc r T , where the scaling exponent r 
is unlikely to depend on the details of the energy supply 
[1]. Below the inertial range, there is a dissipation range. 
The power-law scaling does not persist because of the ex¬ 
istence of the Kolmogorov length rjK as a characteristic 
scale [ 2 ]. 

Turbulence is intermittent at these length scales in the 
inertial and dissipation ranges. An example is the rate of 
energy dissipation per unit mass £ = vY2i j=i (dvt/dxj + 
dvj/dxi) 2 / 2 , where 1 7 is the velocity in the direction of 
the coordinate Xi. This rate is uniform according to the 
1941 theory of Kolmogorov [1], but actually it fluctuates 
in space. If smoothed over a scale in the inertial or dis¬ 
sipation range (see Sec. II), it is significant only within a 
fraction of the space. The fraction decreases with a de¬ 
crease in the scale [3, 4]. Such small-scale intermittency 
has been intensively studied for decades [5]. 

The central issue is to formulate a statistical framework 
that describes the rate e of energy dissipation. Although 
the dissipation itself is negligible in the inertial range, a 
power-law scaling is expected for the rate e r smoothed 
over length scale r in this range, 

(e™)ocr r "\ (la) 

Since the dissipation rate is positive as s > 0, the expo¬ 
nent t m has to satisfy an inequality explained in Sec. II 


for the dimension D = 1, 2, or 3 of the smoothing region 

[ 6 ], 

— > — D at to > 0. (lb) 

TO 

This does not apply to the rate of energy transfer through 
the scale r, which takes both positive and negative values 
[5]. Also from an assumption for the intermittency (see 
Sec. II), we expect an asymptote [7], 

— — > —D as to —> + 00 . (lc) 

TO 

The framework has to be justified by mathematical con¬ 
sistency with these and other constraints [8]. We would 
be left with free parameters that are to be determined by 
physics outside the framework. 

Kolmogorov in 1962 [9] was the first to propose a the¬ 
ory for the small-scale intermittency, by assuming that 
e r obeys a log-normal distribution or equivalently lney 
obeys a Gaussian (normal) distribution throughout scales 
r in the inertial range. However, since the resultant expo¬ 
nent r m does not satisfy Eq. (lb), this log-normal theory 
does not deserve to be a statistical framework of the dis¬ 
sipation rate. 

The log-normal theory is nevertheless of interest be¬ 
cause the probability density function (PDF) of lney is 
self-similar among those scales r. While its position and 
its width are different, its shape remains identical. To see 
mathematics of this self-similarity, we ignore fluctuations 
at the largest scale R of the inertial range and divide the 
scales as r = ri < r 2 < • • • < ryv+i = R, 



The summands ln(£y n /ey n+1 ) are set to be Gaussian and 
to be independent of one another. Since the Gaussianity 
is invariant under addition of its independent variables, 
the sum In (Et/er) is again Gaussian. The self-similarity 
in such a sense, if actual, should be useful to formulating 
the framework. 

Gaussian distributions make up a two-parameter fam¬ 
ily, which is included in a four-parameter family of stable 
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distributions [10-12]. If and only if the distribution is sta¬ 
ble, the shape of its PDF is invariant under addition of its 
variables. For scales r in the inertial range, Kida [13] con¬ 
sidered that In (£ t /£r) obeys a stable distribution. Two 
of its four parameters were determined mathematically, 
while the other two were left as free parameters. This log- 
stable theory is necessary and sufficient for a self-similar 
PDF to describe the fluctuations of In (e r /£R). However, 
the theory as a whole is not consistent with Eq. (lb). In 
fact, its special case is the log-normal theory. 

The same self-similarity of ln(e r /er) exists in the 1941 
theory of Kolmogorov [1], i.e., r m = 0, which is not inter¬ 
mittent but is consistent with Eq. (lb) and is regarded 
as the degenerate case of some stable distribution. We 
search for an extension of this classical theory within the 
whole family of the stable distributions. The result is a 
one-parameter framework where In (st/sr) in the inertial 
range has a self-similar PDF and is consistent with Eq. 
(lb) and also with Eq. (lc). We discuss the application 
to the velocity field, consistency with the actual turbu¬ 
lence, physics that determines the parameter value, and 
differences from the existing frameworks. 


II. BASIC SETTINGS 


The dissipation rate e is smoothed over length scales r 
in the inertial range. Since its boundaries are not definite, 
it is set to be wide enough as compared with ranges of 
scales that could be used as its largest and smallest scales. 
The largest scale is fixed at R, which is well below the 
scale of the energy supply [1]. We do not fix the smallest 
scale because the inertial range could extend to any small 
scale r > r)K in accordance with the Reynolds number of 
the turbulence. 

We consider the one-dimensional smoothing case D = 
1. As illustrated in Fig. 1, the rate e is averaged over a 
segment of length r centered at a position x along a line 
in the three-dimensional space, 

e r {x ) = - f e{x')dx'. (3a) 

X J\x'—x\<r/2 


The direction of this line is fixed arbitrarily in the isotrop¬ 
ic turbulence. At the same position x, we define a random 
variable Xr as 


Xr{x) = In 


re r (x) 

Re R {x) 


< 0 


at r < R. 


(3b) 


The inequality Xr < 0 is due to the positivity e > 0 in 
Eq. (3a). We also consider the three-dimensional smooth¬ 
ing case D = 3, by averaging the rate £ over a spherical 
volume of diameter r centered at a position x, 


and 


£ r (x) = 


~ f 

7r? ’ 3 J\X'-X\<r/2 


£(x')dx', 


(4a) 


Xr(x) = In 


r 3 £ r (x) 

R 3 £r(x) 


< o 


at r < R. 


(4b) 


D= 1 






FIG. 1. Schematic illustration of smoothing regions for D = 1 
in Eq. (3a) and for D = 3 in Eq. (4a) on a two-dimensional 
cut of the space. The grey and dark grey areas denote regions 
of intense dissipation and of very intense dissipation. 


The center of the smoothing region is fixed at a position 
in the space so as to regard £ r and as random functions 
of r in the inertial range up to the largest scale R. Their 
spatial averages are obtained by shifting the position of 
that center. 

The smaller smoothing region is included in the larger 
so long as their centers are common. Also since the rate £ 
is positive, its integral over the region oc r D £ r is increas¬ 
ing in the scale r. The moment ((r D £ r ) m ) oc r Dm+Tm at 
to > 0 is increasing as well, i.e., Dm + r m >0 [6]. This 
corresponds to Eq. (lb), which always holds regardless of 
relative positions of the smoothing regions because the 
turbulence is homogeneous. 

Another constraint is expected from the intermittency 
[7]. We assume that r D £ r /R D £R is independent of R D £r 
and takes any value from 0 to 1. Then r D £ r and R D £r 
take the same maximum value, i.e., ({r D e: r ) m ) 1 / m —> 
{{R D £R) m ) 1 t rn as to —> + 00 . It means r ( Dm + T m.)/ m r o 
and hence Eq. (lc). The maximum values are the same 
even at r -C R. That is, any very intense dissipation has 
been assumed to occur in an isolated point-like region 
(see Fig. 1). Although the energy dissipation is intense 
in organized structures, the intensity within each of the 
structures is not uniform but is enhanced in its point-like 
regions [5]. They are at r > 0, so that they could exist 
in smoothing regions of any dimension D and hence Eq. 
(lc) holds for any value of D < 3. 


III. FORMULATION 

The following conditions are imposed on the random 
variable Xr in the inertial range at r < R: (i) Xr x ~Xr 2 has 
a PDF that depends only on r\jri for any pair of rq and 
Ub (ii) Xri-Xr 2 , Xr 2 ~Xr 3 , XVat-i ~Xr N do not depend 
on one another for any finite series of rq < rq < • • • < r/v; 
and (iii) Xr has a self-similar PDF. If the sign = is used 
to denote that the two random variables have the same 
PDF, the condition (iii) is such that a constant C rii r 2 > 0 
is found for any pair of rq and r2 to have C' ri>r 2 Xr 1 — Xr 2 - 
From these three conditions, it follows that Xr is strictly 
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stable at each of the scales r. The above self-similarity 
is described as 

Xr — [in (.R^/r^)] x* with 0 < a < 2 and /3 > 0. 

( 5 ) 

Here x* is identical to Xr at In (RP/rP) = 1. We are also 
to use an auxiliary condition that er is uniform in space, 
i.e., er = (e), although it is only for simplicity because 
Er does not affect the PDF of Xr nor of e t /er at least 
under the condition (i). 

Here is a mathematical explanation [10, 12]. The con¬ 
ditions (i) and (ii) permit us to regard Xr as a stochastic 
Levy process for the time parameter t = \n(RP/r@) > 0. 
Any process Xt at t > 0 is a Levy process if Xt — 0 at 
t = 0, if Xt is a continuous function of t, if Xti — Xt 2 has a 
PDF that depends only on t\ —t 2 for any pair of t\ and t- 2 , 
and if Xt 1 ~Xt 2 , Xt 2 ~Xt 3 , Xt N = ^Xt N do not depend 
on one another for any finite series of t\ < t 2 < • • • < 

The parameter /3 is to define the scale r that corresponds 
to the time t = 1 . 

The Levy process Xt is stable if constants C tl t2 > 0 
and C'l i t2 are found for any pair of t\ and t 2 to have 
C tl ,t 2 X tl + = Xt 2 • If all the pairs yield C( ut2 = 0, 

the stability is said to be strict. This is our condition (iii). 
The constant Cf t2 serves as a shift of the PDF. It is not 
required for Xt used here, judging from our discussion in 
Sec. II. 

By defining x* as Xt at t = 1, the strictly stable process 
is described as xt — t l ^ a X* with a parameter 0 < a < 2 . 
This is Eq. (5). At each t , the distribution of xt is strictly 
stable with the same parameter a. To be specific, if Xt l 
and Xt 2 are independent copies of Xt> we have C\ x[ \ + 

C 2 Xt ,2 — (Cf+C^) 1/,a Xi for any pair of constants C 1 > 0 
and C 2 > 0. The shape of the PDF is invariant without 
its shift. 

The condition (i) means homogeneity of ln(£ rn /£y n+1 ) 
throughout scales r in the inertial range, while the condi¬ 
tion (ii) means their independence. Under these two, no 
characteristic scale is permissible, so that we expect the 
existence of the power-law scaling (e™) oc r Tm (see Sec. 
I). They are actually idealized conditions. For example, 
ln(e ri /e r2 ) and ln(£ r 2 /e r3 ) are not independent of each 
other for narrowly separated scales, say, r^/ri < 10 [14], 
although such narrow separations are not important if 
the inertial range is wide enough. 

The condition (iii) means a self-similarity of the PDF of 
In {e t /er), albeit requiring its shift, among scales r in the 
inertial range. If the shape of such a self-similar PDF is 
used to define the intermittency, it is not dependent on r. 
The dependence does exist for the intermittency of e t /er 
because the PDF is no longer self-similar. From Eq. (5) 
via Eq. (3b) or (4b), we find that E r j er is increasingly in¬ 
termittent with a decrease in the scale r from the largest 
scale R. 

Stable distributions are described generally by four pa¬ 
rameters, but those in the strict sense make up a family of 
three parameters, where excluded is a parameter to shift 
the PDF [10, 12]. Except for a — 1, the characteristic 


function of the distribution of x* is 

<exp(ix.O) = exp (-A|ere^/ 2 I«I) . (6) 

While A > 0 determines the width of the PDF, a and 9 
determine its shape (|#| < a for 0 < a < 1 and |f?| < 2 — a 
for 1 < a < 2). The stable distributions for a = 2 and 
9 = 0 are Gaussian. 

We focus on the distributions for 0 < a < 1 and 9 = a. 
They alone are totally skewed to the left, i.e., x* < 0 
[10, 12], which is required from Eq. (3b) or (4b) via Eq. 
(5). By using the above values of a and 9 in the inverse 
Fourier transform of Eq. (6), the PDF /(x*) is written 
as 


0 


at x* > 0 


fix*) 


E 


A(—A)" 1 T(na + 1) sin( 7 rna) 

7 m !| x *| n “ +1 


at x* < 0 . 


( 7 ) 

Here T is the gamma function. The tail for x* ~^► — oo 
is in the form of power law oc |x*| - ^“ +1 \ which leads to 
divergence of moments such as the average (x*)- 

To formulate the one-dimensional smoothing case D = 
1, we use Eq. (3b) to have (r m s™ / R m e 1 ^) = (exp (mx r ))- 
Then (exp(?nXr)) is obtained by substituting Eq. (5) into 
Eq. ( 6 ) and by replacing £ with — im. This replacement 
corresponds to an analytical continuation from £ to —wr¬ 
it is justified at in > 0 by the convergence of (exp(mx*)), 
which is in turn justified by the form of fix*) i n Eq. (7). 
The result is 



(r/ R) ~ m + rn °‘^ x exp(7ra/2) at m > 0. (8) 


However, does not exist at m < 0. The reason is 

the power-law tail of f(x*) and the resultant divergence 
of (exp(rrax*)) at to < 0. Negative-order moments ob¬ 
served in numerical simulations and in experiments are 
attributable to, e.g., random noise that is inherent in any 
actual data. 

Finally, the parameters (3 and A are eliminated by mak¬ 
ing use of the auxiliary condition er = (e) . The left-hand 
side of Eq. ( 8 ) is reduced to (£( ra )/(£) m . At m = 1, it is 
(Er) / (e) = 1. The corresponding exponent for r/R in the 
right-hand side is —1 + /3A exp( 7 ra/ 2 ) = 0. Accordingly, 


(fn 

(e) m 


(r/R)- m+m “ , 


and hence as the exponent r TO for D = 1, 


(9a) 


T. m = —in + m a with 0 < a < 1 at to > 0. (9b) 

Even if er is not uniform, er does not affect e t /er under 
the condition (i). We have {e™/£r}{£r) = (e™), through 
which Eq. ( 8 ) leads to (e™) oc r Tm with the same Eq. (9b). 
Then, Eq. (9a) could be regarded as a result conditioned 
on a given value (e) of the fluctuating rate Er [4]. 

This is the necessary and sufficient framework for the 
case of D = 1 where In {Et/er) has a self-similar PDF 
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and satisfies Eq. (lb). Also satisfied is Eq. (lc). The in- 
termittency is described by the parameter a. Especially 
in the limit a —> 1, we reproduce the 1941 theory of Kol¬ 
mogorov [1], i.e., r m = 0, for which Eq. (lc) does not 
hold because the intermittency does not persist. 

The three-dimensional smoothing case D = 3 is also 
formulated, by replacing the length ratio r/R in Eq. (9a) 
with the volume ratio r 3 /R 3 , 

= (r*/R 3 ) - m+m “ = (r/R)~ 3m+3ma , (10a) 

and hence as the exponent r m for D = 3, 

r m = —3m + 3 m a with 0 < a < 1 at m > 0. (10b) 

Between D = 1 and D = 3, the exponent r m is different. 
We have r m = —Dm + Dm a , which holds for D = 2 as 
well. The value of the parameter a could depend on the 
smoothing dimension D, although there has been found 
no mathematical relation. 

The present formulation is almost equivalent to that of 
Kida [13]. However, his stable distribution for \ r is not 
in the present three-parameter strict sense but is in the 
four-parameter general sense. The resultant exponent 
has the form of T m = —fi(m a — m)/(2“ — 2) with two free 
parameters fi > 0 and 0 < a < 2, even if the smoothing 
dimension D is given. We have imposed an additional 
constraint of Xr < 0 to eliminate the parameter /i and to 
obtain Eqs. (9b) and (10b) where r TO always satisfies the 
inequality of Eq. (lb). 


IV. APPLICATION TO VELOCITY FIELD 

The power-law scaling in the inertial range is expected 
also for the two-point velocity difference Su r (x ) = u{x + 
r) - u(x), 

(|(b4 ra |) a at m > 0. (Ha) 

Here u is longitudinal or lateral, i.e., parallel or perpen¬ 
dicular to the line through the two points. The coordi¬ 
nate x is along this line. We study instead of Su r 
itself [9]. Since 5u r could take the value of 0, the moment 
(|<5u™|) does not exist at m < 0. 

If the velocity u is to be bounded such that the flow 
does not become supersonic and hence does not become 
compressible, there is a necessary condition 

^>0. (lib) 

dm 

This is because l\8u™ +m ' |)/(|<5u™|) oc has to 

be bounded for m! > 0 even if the inertial range extends 
to any small scale, i.e., £m+m' ~ Cm > 0 [15]. 

The exponent is obtained from the exponent r m for 
the dissipation rate e r . Based on a velocity scale ( re r ) 1 ' 3 , 


there is a relation known as Kolmogorov’s refined simi¬ 
larity hypothesis [9], 

{\6v?\) oc ((re r ) m/3 ) and Cm = T m/3 + y. (12) 

To the longitudinal velocity, Eq. (12) is plausible because 
the energy dissipation is intense in regions under intense 
strain. We do not use the lateral velocity. Its two-point 
difference is rotational and is not related closely to the 
rate of the dissipation [16]. 

We consider the one-dimensional smoothing case D = 
1. The two points for the velocity difference 5u r are the 
ends of the line segment used to smooth the dissipation 
rate e r . From r TO = —m + m a in Eq. (9b) via Eq. (12), 

( Tfi \ ^ 

—) with 0 < a < 1 at to > 0. (13) 

The exponent Cm satisfies Eq. (lib) and other constraints 
such as C 2 > 2/3 [17]. In the limit a —» 1, we reproduce 
the 1941 theory of Kolmogorov [1], i.e., Cm = m/3. 

We also consider the three-dimensional smoothing case 
D = 3. The two points for the velocity difference Su r are 
poles of the sphere used to smooth the dissipation rate e r . 
From T m — — 3m+3m“ in Eq. (10b) via Eq. (12), we have 
Cm = — 2m/3 + 3(m/3)“. As m —> + 00 , the exponent Cm 
becomes inconsistent with Eq. (lib) and becomes nega¬ 
tive. The reason is r m /m —> —3 that satisfies Eq. (lc). 
Thus, Eqs. (lc), (lib), and (12) are inconsistent with one 
another in the case of D = 3. 

The dissipation exponent r m is different between D = 
1 and D = 3, while the velocity exponent Cm is unique. It 
is hence possible that the relation of r m to Cm in Eq. (12) 
is suited to D = 1 but is not to D = 3. The velocity dif¬ 
ference Su r is an integral of the velocity derivative du/dx 
over the separation r, while an integral of [du/dx) 2 is one 
of the components of the dissipation rate e r . Regions of 
these two integrations are identical in the case of D = 1. 
They are not identical in the case of D = 3, for which it 
would be of interest to reconsider whether or not any re¬ 
lation exists between the velocity difference Su r and the 
dissipation rate s r . 

V. COMPARISON WITH OBSERVATIONS 

Figure 2 compares the scaling exponents of our frame¬ 
work with those observed in numerical simulations of ho¬ 
mogeneous and isotropic turbulence [18-20]. Since our 
framework is well consistent with the simulations, we con¬ 
clude that our conditions (i)-(iii) used in Sec. Ill are at 
least good approximations of the actual turbulence. This 
conclusion applies especially to the condition (iii), which 
says that ln(e r /£R) has a self-similar PDF. 

The three-dimensional smoothing case D = 3 has been 
studied for the dissipation exponent r m in Eq. (10b). It 
yields a ~ 0.90-0.95. On the other hand, in the one¬ 
dimensional smoothing case D = 1, the data of r m are not 
available, except for a surrogate of the dissipation rate 
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FIG. 2. Comparison of our framework with numerical simu¬ 
lations of forced steady states of homogeneous and isotropic 
turbulence at the Reynolds number for the Taylor microscale 
of 216 [18], 460 [19], and 600 [20]: (a) dissipation exponent 
Tm for D = 3 in Eq. (10b), and (b) velocity exponent for 
D = 1 in Eq. (13). The dotted lines denote the 1941 theory 
of Kolmogorov [1]. 


oc ( du/dx ) 2 that tends to be more intermittent [21]. We 
have instead used the velocity exponent ( m in Eq. (13). It 
yields a ~ 0.7-0.8. Thus a takes different values between 
the cases of D = 1 and D = 3. 

With a numerical simulation or with an experiment, 
we would have to compare our framework again in the 
future. The conditions (i) and (ii), i.e., homogeneity and 
independence of ln(e rn /e rn+1 ) in the inertial range, would 
hold if the Reynolds number were high enough to set up 
a wide inertial range. This is not the case in the existing 
simulations and experiments [5, 22]. We expect that the 
future data are at a high enough Reynolds number, which 
would also permit the more accurate observations of r m 
and ( m . The auxiliary condition, i.e., uniformity of £r. 
also does not hold because £r fluctuates significantly even 
if the turbulence is fully developed and is filling the space 
[23, 24]. Although such fluctuations are unlikely to affect 
the power law in the form of (e™) oc r Tm , we expect that 
the future data are of enough size to study the scaling of 
(e™) at each value of £r. Then, our framework should 
exhibit a better consistency with the simulation or with 
the experiment if our self-similarity condition (iii) holds 
for the actual turbulence. 


VI. EXPLANATION FOR PARAMETER VALUE 

Thus far, we have discussed mathematics of the dissi¬ 
pation rate e r to formulate a statistical framework of its 
log-stable law. If this framework does apply to the small- 
scale intermittency, the observed value of the parameter 
a is explainable by its physics. Here is an attempt of 
such an explanation, which is based on the transition 
from the inertial range to the dissipation range in the 
one-dimensional smoothing case D = 1. 

The transition is determined by the Reynolds number 
Re r . For a region of length r centered at some position 
x , we follow Kolmogorov [9] to define Re r by making use 
of a velocity scale ( re r ) 1 / 3 , 


Re r (a;) 


r[r£ r (a:)] 1 / 3 

v 


(14a) 


The same velocity scale exists in Eq. (12). From Re^ = 1, 
the transition length scale rj is defined as 


V(x) 


t/3/4 

4 /4 0) 


(14b) 


The scale ry fluctuates spatially with fluctuations of the 
dissipation rate e v [9]. Since (ry _4m ) is related to (e™), we 
assume that (e™) in the inertial range transits to (e™) in 
the dissipation range at the length of r = C'*(? 7 _4m )~ 1 / 4m 
[25]. The constant (7* > 10 is due to the Reynolds num¬ 
ber Re r = C*/ 3 for the actual occurrence of the transition 
[5], 

The inertial range exhibits (e™)/{e) m = ( r/R) Tm if we 
use Eq. (9a), which yields 



?=ln(r /R). 


(15a) 


Here dT m /dm = a/m 1-a — 1 > 0 at to < to* = a 1 /! 1- ") 
and dT m /dm < 0 at to > to*. The value of to* is found 
at 0.3-0.4 in the case of a > 0.5. 

The dissipation range exhibits smooth velocities such 
that dvi/dxj ~ constant in any region of length r < (7*77, 
from which we have e r (x) ~ £ v (x) and (s™) ~ {e™). Also 
by using the Kolmogorov length t/k = v 4 I 4 /{e) 4 I 4 [1], 




TO 




TO 



n-’b* 

The derivative has been approximated by a linear slope 
from In ((e 7 / 1 )/(e) m ) = 0 at to. = 0. It holds if to is small, 
say, m < to*. If m were not small, the slope would have 
to be taken from ln((£™)/(s)™) = 0 at m = 1. 

Finally, our assumption is used to equate Eq. (15a) 
with Eq. (15b) at r = C^(r]~ 4m )~ 1 ^ 4m . The result for 
to = 1 /4 is 


R 

C*r] 


\ 4/(4+dr m /dm )\ m=x/i 

<7*?7k ) 


(16) 
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Here ( R/C * 77 ) serves as the width of the inertial range. 
The classical width R/C^ij k is reproduced in the limit 
a —> 1, corresponding to the 1941 theory of Kolmogorov 
[ 11 - 

Given the classical width of the inertial range R/C^rjK, 
which is fixed by conditions external to the turbulence, 
(R/C^rj) is minimal at a = 0.72. Since a ~ 0.7-0.8 is ob¬ 
served for D = 1 in Fig. 2, we consider that the observed 
value of a corresponds to the minimum of the width of 
the inertial range (R/C*rj). Its inverse (D/C*??) -1 serves 
as the mean efficiency of the energy dissipation. While 
is large in regions where the kinetic energy is effectively 
dissipated into heat, is not large in the other regions. 
The average over all the regions is (? 7 -1 ) -1 . Thus, a plau¬ 
sible explanation is that the value of a is determined so 
as to maximize the mean efficiency of the energy dissi¬ 
pation (R/C^rj)- 1 and thereby to minimize the width of 
the inertial range (R/C*rf). 

This explanation has used our auxiliary condition £r = 
(e) and the resultant Eq. (9a). Even if £ R is not uniform, 
the same explanation is possible by considering R/C*r /k 
and (R/C^rj) for each value of £r (see Sec. III). 

To study the three-dimensional smoothing case D = 3, 
we have to consider (R 3 /(C^r/) 3 ) instead of {R/C^rf), i.e., 
m = 3/4 instead of m = 1/4. However, even if Eq. (15b) 
is rewritten for m > m*, it leads to an inconsistent result 
that (R 3 /(C*r]) 3 ) is not minimal in the range of 0 < a < 
1. The reason is attributable to the velocity scale (rey) 1 / 3 
in Eq. (14a), which is unlikely to be suited to D = 3 as 
discussed for Eq. (12) in Sec. IV. An alternative approach 
is required to explain the observed value of the parameter 
a and the relation to its value for D = 1, if possible on 
the basis of some definition of the mean efficiency of the 
energy dissipation. 


VII. COMPARISON WITH OTHER THEORIES 

For the small-scale intermittency, most of the existing 
theories are based on the conditions (i) and (ii) used also 
in Sec. III. They yield (e?/e%) = {e™ R /e r £) XlIl{r/R)/hx * 1 
with r < jR < R, which is necessary and sufficient for 
the existence of a power law (e™/£ R ) = ( r/R) Tm [6]. It is 
equivalent to (el/ 1 ) oc r Tm because £ t /£r does not depend 
on e r under the condition (i). To constrain the exponent 
r m , each theory has some additional condition. 

The current paradigm is multifractality [26]. To for¬ 
mulate this into a statistical framework [27, 28], instead 
of using our self-similarity condition (iii), the large devia¬ 
tion theory is used for the sum of In (e r „Ar„ +1 ) in Eq. (2) 
with r n /r n +i = constant, i.e., with N oc — ln(r/i?). This 
limit theorem says that the PDF of ln(e r /£ij)/In (r/R) 
has a specific form of asymptote in the limit r/R 0 [27 
29]. The inertial range is assumed to extend to that limit, 
where the asymptotic PDF is related to the exponent r m 
via a series of differences between the smoothing dimen¬ 
sion D and some fractal dimensions on the D-dimensional 
cut of the space. Since these differences are independent 


of the value of D [11], the exponent r m is independent as 
well. 

For the multifractal PDF of ln(£ rtl /£ rn+1 ), there have 
been proposed various theories. Some of them lead to 
(£ 7 / l /£ R ) even at m < 0 [30]. The others do not as in our 
case of the stable distribution [31]. 

To any stable distribution that is not Gaussian, the 
large deviation theory is not applicable. It requires that 
the tail of the PDF of ln(e r /e R )/ln{r/R) decays expo¬ 
nentially as ln(r/.R) —> — 00 , i.e., as N — > + 00 . This is 
not the case for the PDF of the stable distribution that 
retains its power-law tail [29]. Especially for 0 < a < 1 
studied here, ln(£ r /e R )/ \a(r/R) diverges in that limit as 
TV 1 /" -1 (see Sec. III). The fractal dimensions are hence 
not available. 

The above discussion means that any multifractal PDF 
of ln(£ rn /£ rn+1 ) is not stable and is not self-similar. An 
example is the Poisson distribution [31 -33]. Its param¬ 
eter is A > 0. If ln(£ ri /£ T . 2 ) and ln(£ r2 /£ r3 ) are related 
to the parameter values A\ and A 2 , their sum \n(£ ri /e T3 ) 
is related to a different parameter value A\ + A 2 , which 
results in a different shape of the PDF. 

The independence of the scaling exponent r m from the 
smoothing dimension D is an important characteristic of 
any multifractal PDF. Among the dimensions D = 1, 2, 
and 3, the exponent r TO is unique. Since it has to satisfy 
Eq. (lb) for D > 1, it satisfies Eq. (lc) at most for D = 1, 
although we consider that r m is not unique and satisfies 
Eq. (lc) for each value of D as in our case of the stable 
distribution. 

Multifractality is thus distinct from our log-stability. 
The underlying assumption for using the large deviation 
theory is that the intermittency is asymptotically deter¬ 
mined in the course of the mean energy transfer to the 
smaller scales, i.e., the energy cascade [34]. Nevertheless, 
the energy is transferred locally and instantaneously to 
the larger scales as well as to the smaller scales [5]. These 
scales interact with one another and should have settled 
into a self-similar state that is to be described by a stable 
distribution of In (e r /£ii)- In addition, it is too idealized 
to ignore the existence of the dissipation range by taking 
the limit r/R —> 0 for scales r in the inertial range. If 
the limit is not taken, the large deviation theory is not of 
use. The exponent r m remains available from the power- 
law scaling itself, but it has to be positioned within some 
alternative to the multifractal framework. 

The same discussion applies to other frameworks such 
as those based on the central limit theorem [35] and on 
the generalized central limit theorem [36], aside from the 
consistency or inconsistency with Eq. (lb) and so on. 

With the exponents r m and observed thus far in nu¬ 
merical simulations and in experiments, the multifractal 
and other theories exhibit a consistency almost compara¬ 
ble to that of our log-stable law [18-20]. The theories are 
not judged definitely at this level of observations. As for 
the observational judgment between the log-stable law 
and the others, it is essential to study whether or not 
the PDF of ln(£ r /£/j) is self-similar among scales r in a 
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wide enough inertial range. This is to be expected in the 
future. 


VIII. CONCLUDING REMARKS 

The small-scale intermittency has been studied about 
a fully developed state of three-dimensional homogeneous 
and isotropic turbulence in an incompressible fluid. For 
the rate £ r of energy dissipation smoothed over length 
scale r in the inertial range up to the largest scale R, 
we have followed Kida [13] to impose a self-similarity on 
the PDF of ln(s r /en). Such a self-similarity is expected 
from interactions among the scales through the local and 
instantaneous transfer of the kinetic energy. The result 
is an extension of the 1941 theory of Kolmogorov [1], i.e., 
a one-parameter framework where \n(e r /£ r) obeys some 
stable distribution. As a statistical framework of the self¬ 
similar PDF of In (et/sr), this log-stable law is necessary 
and sufficient. 

The scaling exponents have been obtained for the scal¬ 
ing of the dissipation rate {e™) oc r Tm in Eqs. (9b) and 
(10b) and for the scaling of the two-point velocity differ¬ 
ence (|<5w™|) oc r < ’ m in Eq. (13). They are consistent with 
theoretical constraints such as Eqs. (lb), (lc), and (lib) 
and with the observed exponents in Fig. 2. 

We have attempted to explain the observed value of 
our parameter a. This is also to explain the existence of 
organized structures in the turbulence because our frame¬ 
work is just a statistical description of them. In fact, the 
observed value of a is distinct from the value a —> 1 for 
the classical case of no structures, i.e., the 1941 theory of 
Kolmogorov [1], By using Eq. (16), we have related the 


observed value of a to the minimum width of the inertial 
range and thereby to the maximum mean efficiency of the 
energy dissipation. It follows that the intermittency oc¬ 
curs and the structures are organized so as to maximize 
the mean efficiency of the dissipation. 

We have not considered the fluctuations of £r at the 
largest scale R of the inertial range. They do not have to 
be considered if the power-law scaling (e™) oc r Tm is ex¬ 
act. It means that the inertial range is not characterized 
by any of the scales, including scales of those fluctuations. 
This is the case of our framework where £ t /sr is set to be 
independent of £r. but we would have to study the very 
case of the actual turbulence. Oboukhov [37] proposed a 
PDF of those fluctuations. Since it has turned out to be 
a good approximation [24, 38], it would be useful to that 
study. 

The stable distribution is also expected for a logarithm 
of any other positive random variable, e.g., enstrophy [5, 
18], if its PDF is self-similar in a range of the scales. This 
is not limited to variables of the turbulence. For example, 
a power-law scaling is found in a variety of fields [11] such 
as the cosmological density field observed through the 
galaxy number density [39]. They have been discussed to 
be fractal or multifractal, but some of them could be log- 
stable. It is of interest to search for the stable distribution 
among logarithms of these and other positive random 
variables. 
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